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ABSTRACT 

The Bethe-Salpeter (BS) equation for scalar-scalar bound states in scalar the- 
ories without derivative coupling is formulated and solved in Minkowski space. 
This is achieved using the perturbation theory integral representation (PTIR), 
which allows these amplitudes to be expressed as integrals over weight functions 
and known singularity structures and hence allows us to convert the BS equa- 
tion into an integral equation involving weight functions. We obtain numerical 
solutions using this formalism for a number of scattering kernels to illustrate 
the generality of the approach. It applies even when the naive Wick rotation 
is invalid. As a check we verify, for example, that this method applied to the 
special case of the massive ladder exchange kernel reproduces the same results 
as are obtained by Wick rotation. 



1. Introduction 



There has been considerable recent interest in covariant descriptions of bound 
states, for example, in conjunction with model calculations of high-energy processes 
such as deep inelastic scattering. A fully covariant description of composite bound 
states is essential for the understanding of hadronic structure over the full range of 
available momentum transfer. In a relativistic field theory the two-body component 
of a bound state is described by the appropriate proper (i.e., one-particle irreducible) 
three-point vertex function or, equivalently, by the Bethe-Salpeter (BS) amplitudes. 
The BS amplitude satisfies an integral equation whose kernel has singularities due to 
the Minkowski metric. The resultant solutions are not functions but mathematical 
distributions. The singularity structure of these distributions makes it difficult to 
handle the BS equation numerically in Minkowski space. 



One approach to dealing with the difficulties presented by the Minkowski-space 
metric is to perform an analytic continuation in the relative-energy variable p°, which 
is the so-called "Wick rotation" @. This has been widely used as a means of solving 
model BS equations, e.g., see Ref. 3 and references therein. The special case of the 
ladder BS equation is solved as a function of Euclidean relative momentum in the 
standard treatment. However it is necessary to rotate Euclidean solutions back to 
Minkowsi ones in order to apply the BS amplitude (vertex function) to calculations 
of absorptive part of scattering amplitudes for example. The difficulties associated 
with this approach in the general case arise from the fact that since the total four- 
momentum P must remain timelike then the kernel becomes complex, so that one 
needs to perform analytic continuation numerically to obtain Minkowski solutions. 
In addition, when one uses a "dressed" propagator for the constituent particles or 
more complicated kernels in the BS equation, the validity of the Wick rotation be- 
comes highly nontrivial. For example, essentially all of the dressed propagators stud- 
ied previously in Wick-rotated Dyson-Schwinger equation approach contain complex 
"ghost" poles!. Hence, the naive Wick rotation obtained by the simple transcription 
of Minkowski metric to the Euclidean metric and vice versa is not valid in general. 
So while the Euclidean-based approaches certainly play a very important role and 
can be useful in model calculations, it is preferable to formulate and solve the BS 
equation directly in Minkowski space. Here we present a method to solve the BS 
equation without Wick rotation by making use of the perturbation theory integral 
representation (PTIR) for the BS amplitudeS. 

The PTIR is a natural extension of the spectral representation for a two-point 
Green's function to an n-point function in a relativistic field theory. Since the Feyn- 
man parametric integral always exists for any Feynman diagram in perturbative cal- 
culations, one can always define the integral representation such that the number of 
independent integration parameters is equal to that of invariant squares of external 
momenta. Each Feynman diagram contributes to the weight distribution of the para- 
metric integral, so that the complete weight function for the renormalized n-point 
function is identical to the infinite sum of Feynman diagrams for the renormalized 
Lagrangian of the theory. Hence, we see that the PTIR of a particular renormalized 
n-point function is an integral representation of the corresponding infinite sum of 
Feynman diagrams for the renormalized theory with n fixed external lines. 

The PTIR method was first used in conjunction with a Wick rotation by Wick 
and Cutkosky for solving a scalar-scalar bound state with a massless scalar exchange 
in the ladder approximation!'!. This is now commonly referred to as the Wick- 
Cutkosky model. They solved the BS equation in terms of a single variable integral 
representation, which is a special case of the PTIR. After initial efforts to solve scalar- 
scalar BSE's using, the PTIR methodS, Nakanishi made a detailed and systematic 
study of the PTIRb. Here we extend the PTIR method to solve the BS equation with 
a broad class of scattering kernels, by transforming the BSE in momentum space to 



an integral equation for the weight function without relying on the Wick rotation. 



2. Scalar-Scalar BS Equation 

Let us consider a bound state of two spinless particles having a mass m. They 
interact with each other through the exchange of other spinless particles. The Bethe- 
Salpeter amplitude P) for the bound state having the total momentum P and 
the relative momentum p obeys the following equation: 

[D(P/2+p)D(P/2-p)]- 1 ^p,P) = J -Al-/(p,g;P)$(g,P) (1) 

where D(q) is the propagator of 0-particle. We approximate it with the free one: 
D°(q) = 1/ (m 2 — q 2 — ie). The scattering kernel I(p, q; P) describes the process <fxp — > 
<fxft, where p and q are the initial and final relative momenta respectively. It is given by 
the infinite series of Feynman diagrams which are two-particle irreducible with respect 
to the initial and final pairs of constituent <fi particles. For purely scalar theories 
without derivative coupling we have the formal expression for the full renormalized 
scattering kernel 

i(p,g;P) = l<hld£\ r ^ 4 ^4 . : , . (2) 
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where q 2 is the 4-momentum squared carried by fa and s,t and u are the usual 
Mandelstam variables. This expression has been derived by Nakanishi (PTIR)i. Since 
only six of these are independent due to the relation qf + q^ + qi + ql = s + t + u, this 
is also the number of independent ^-parameters. Hence, we need only introduce the 
"mass" parameter 7 and six dimensionless Feynman parameters £j with a constraint. 
The symbol Q denotes the integral region of & such that ft = | < & < 1, Y^d — 
l(i = 1, . . . ,6)}. The "mass" parameter 7 represents a spectrum of the scattering 
kernel. The function p c h{l,0 gives the weight of the spectrum arising from three 
different channels which can be denoted ch = {st}, {tu}, {us}. This PTIR expression 
follows since any perturbative Feynman diagram for the scattering kernel can be 
written in this form, and hence this must also be true of their sum. In a perturbative 
calculation the weight function p ch is calculable as a power series of the coupling 
constant for a given interaction Lagrangian. 



3. PTIR for BS Amplitude 



Let us now consider the PTIR of the BS amplitude itself. Since the BS amplitude 
in the center-of- momentum rest frame [i.e., with P M = (M, 0)] forms an irreducible 
representation of the 0(3) rotation group, we can label all of the bound states with 
the usual 3-dimensional angular momentum quantum numbers i and l z . We can thus 
construct the integral representation of the partial wave BS amplitude in the rest 
frame and the PTIR automatically allows us to then boost to an arbitrary frame. 

For simplicity, let us begin by considering the s-wave amplitude. From general 
analysis of the PTIR for a one-particle-irreducible 3-point vertex function^, one can 
write the BS amplitude for a s-wave bound state 

oo 1 

Hp, P) = -i [ da I dz— _ ( (3) 
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m 2 + a — [p 2 + zp ■ P + — ie 



where the "mass" parameter a represents a spectrum of the BS amplitude and the 
weight function (p n (a,z) gives its weight. One may regard the dependence of the 
weight function on the dimensionless parameter z as kinematical effects. This two- 
variable PTIR is valid for the vertex function with a timelike total momentum satis- 
fying < P 2 < Am 2 . We have introduced a non- negative integer dummy parameter 
n. We see that a partial integration of Eq. (§) with respect to a connects weight 
functions with different n, i.e., 



y? n +i(a, z) = (n + 2) / da'(p n (a',z) , (4) 

— oo 

provided that we have the boundary conditions 

(fi n (a = —oo,z) = and lim ^ n ^ a ' I = Q (5) 

which must be satisfied in order for the expression (|]) to be meaningful. While 
the non-negative integer n is arbitrary, a judicious choice can be advantageous in 
numerical calculations, since the larger we take n, the smoother the weight function 
becomes for a given P). 

The extension of the PTIR to the BS amplitudes for higher partial waves, i.e, 
those for bound states with non-zero angular momentum [I > 0) is straightforward. 
The momentum-dependent structure, i.e., the denominator, of the PTIR's for the 
bound state BS amplitude, is independent of transformations under the little group 
belonging to the 3- momentum p in the rest frame of the bound state, i.e., P M = (M, 0). 
Hence, for higher partial waves with angular momentum quantum number I and 
third component £ z , $ can be written as the product of the £-th order solid harmonic 
y^ip) = \p\ e Ye z (Pl an d the corresponding PTIR for the scalar (i.e., s-wave) bound 
state in this frame 0. By considering the Lorentz boost from the rest frame to an 



arbitrary one we have the integral representation of the partial wave BS amplitude 
*M(p,P) = -iy e e°(p>) [ da [ dz— „ — , ]W+2 , (6) 
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[m 2 + a - (p 2 + zp ■ P + P 2 /4) - ie] 



where P is an arbitrary timelike 4- vector with P 2 = M 2 and p' = A~ 1 (P)p. The 
Lorentz transformation A(P) connects P and the bound-state rest frame 4- vector 
P' = (M, 0), i.e, P = A(P)P'. In the following sections we will study the BS equation 
Eq. (JJ) in an arbitrary frame in terms of this integral representation. 

4. BS Equation for the Weight Function 

In this section we will reformulate the BS equation Eq. ([!]) as an integral equation 
in terms of the weight functions. This is the central result of our work. For notational 
convenience we first rewrite the scattering kernel I(p,q; P) as follows 

i(p,<r,P) (7) 

' Pch(7>Q 
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where {a c h, b c ^, c c h, . . . , / c h} are different linear combinations of £j in each of the three 
channels, ch— {st}, {tu}, {us}. Substituting the integral representation of the par- 
tial wave BS amplitude Eq. @, together with the PTIR for the scattering ker- 
nel Eq. (H), into the BS equation Eq. ([!]) and after multiplying the propagators 
DiPjl + p)D(P/2 — p) into both sides, one can perform the g-integral without rely- 
ing on Wick rotation imposing the condition for the dummy parameter n such that 
n + 1 > 1/2. This ensures that the g-integral is finite. Using the uniqueness theorem 
of PTIPJ, we obtain the following integral equation for ipw(a, z): 

00 1 

tp$(a,z) = / da / dz)C$(a,z;a,z)(p$(a,z) . (8) 
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We have defined the total kernel function fC^(a, z; a, z), which is a superposition of 
the kernel functions K$(a, z\ a, z), which in turn are defined with a fixed parameter 
set of the scattering kernel {a c h, 6 c h, c c h> • • • , /ch}, 

K$(a,z',a,z) = Y,J d ^J rf fPch(7, D K n(®, *5 «> z )- ( 9 ) 
ch a 

The kernel function with a fixed kernel parameter set is 

*«(«,*«„) = ^\ (-D'jf • wS^S'Zcr 1 (10) 
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with functions A(a,z), B(a,z), R(z,z) and G(z;y) as 



p2 

A(a, z) = a + m 2 — (1 — z2 )—^~i 

( P 2 \ P 2 
B(a, z) = aa + 7 + (a — c) ( m 2 — j — (Ad — 2ez) — , 

q/ + A 

^ = \^6(z -z) + \^6(z - z), 
v ' ' 1 - z K ' 1 + z v y 

and the constants; C = 07 — A (m 2 — — (Aad— e 2 )^-, A = ac—^j. The integration 
over the Feynman parameter y should be carried out using the Hadamard finite part 
prescription: 

Pf ■ — = limRe- — , for integer n. (12) 

x n <^o (x ± ie) n 

This is consistent with the ordinary ie prescription for a calculation of Feynman 
diagrams in momentum space § and we keep e finite but small for numerical studies 
to regularize possible singularities of the kernel. A detailed analysis of singularities 
of the kernel function may be found in Ref. 6. 

Since the weight functions p c h.(j, Q for the scattering kernel are real functions by 
their construction, the total kernel function K$(a, z\ a, z) is real, so that the Eq. (||) 
is a real integral equation with two variables a and z. Thus we have transformed 
the BS equation, which is a singular integral equation of complex distributions, into 
a real integral equation which is frame- independent. As is evident from Eq. (§) if n 
is chosen sufficiently large these can be always transformed into arbitrarily smooth 
functions suitable for numerical studies. 

5. Numerical Studies 

To illustrate the approach we have numerically solved the following two explicit 
examples : 

(a) Scalar-scalar ladder model with massive scalar exchange: The simple ^-channel 
one-a-exchange kernel is given by 

I(p,q;P)= 2 / 2 . (13) 

m a — (p — q) t£ 

The BS equation with this kernel together with perturbative propagator D° is often 
referred to as the "scalar-scalar ladder modef'S. The Wick rotated BS equation for 
this kernel has been studied numerically^. We use this kernel as a check of our 
calculations. 



(b) Generalized ladder kernel: A sum of the one-cx-exchange kernel Eq. (|13| ) and a 
generalized kernel terms in the st-channel with constant kernel parameter sets 
weighted with a factor (l/4)[g 2 /(4ir) 2 ]. We generated {£, st } by a random number 
generator. After the Wick rotation this kernel becomes complex due to the p ■ P 
and q ■ P terms, so that solving the BS amplitude as a function of Euclidean relative 
momentum would be difficult in this case. 

The scattering kernel (a), i.e., the one-cr-exchange kernel corresponds to choosing 
for the kernel in Eq. (|3|) say: p tu = p us = and in the st-channel 7 = m 2 , and 
a st — c s t — 1) b st = —2, d st = e st — f s t — with p st to be some appropriate product of 
(^-functions multiplied by g 2 . It is often convenient (and traditional) to factorize out 
the coupling constant g 2 together with a factor of (47r) 2 , and define the "eigenvalue" 
A = g 2 1 (47r) 2 . Then we solve the equation as a eigenvalue problem with a fixed bound 
state mass P 2 . To regularize singularities of kernel function we kept the regularization 
parameter e in the Hadamard finite part Pf-l/x n small but finite, typically ~ 10~ 5 or 
less. We began by discretizing the a and z axes and then solved the integral equation 
by iteration from some initial assumed weight function. We verified that our solutions 
were independent of e provided it was chosen sufficiently small. It was also, of course, 
confirmed that the solution was robustly independent of the choice of starting weight 
function, the number of grid points, and the maximum grid value of a (when the 
latter two of these were chosen suitably large). We found that the convergence and 
the stability of the eigenvalues with varying the number of grid points, e, and a max are 
excellent. We have solved over a range of P 2 and m a and all solutions reproduced the 
eigenvalues obtained in Euclidean space after a Wick rotation by Linden and Mitter 
@. For a nontrivial example of the scattering kernel, we studied the model kernel (b) 
symmetrized in z — > —z so that we obtain appropriate symmetry for s-wave normal 
solutions. We found that the convergence and the stability of the eigenvalues with 
varying the number of grid points, e, and a max are excellent as for the pure ladder 
case i. Additional details and figures can be found elsewhere i. 

6. Summary 

We have derived a real integral equation for the weight function of the scalar- 
scalar Bethe-Salpeter (BS) amplitude from the BS equation for scalar theories with- 
out derivative coupling. This was achieved using the perturbation theory integral 
representation (PTIR), which is an extension of the spectral representation for two- 
point Green's functions, for both the scattering kernel [Eq. @] and the BS amplitude 
itself [Eqs. (§)]. The uniqueness theorem of the PTIR and the appropriate applica- 
tion of Feynman parameterization then led to the central result of the paper given 
in Eq. (§). We demonstrated that this integral equation is numerically tractable for 
both the pure ladder case and an arbitrary generalization of this. We have verified 
that our numerical solutions agree with those previously obtained from a Euclidean 



treatment of the pure ladder limit. 

This represents a potentially powerful new approach to obtaining solutions of the 
BS equation and additional results and applications are currently being investigated. 
The separable kernel case should be studied, since it has exact solutions and so is a fur- 
ther independent test of the approach developed here. While our detailed discussions 
were limited to the equal mass case, it is worthwhile to investigate generalizations 
including nonperturbative constituent propagators and the heavy-light bound state 
limit to see under which conditions an approximate Klein-Gordon equation can result. 
It is also important to find a means to generalize this approach to include derivative 
couplings and fermions. 
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